(12) NACH DEM VERTRAG UBER DIE INTERNATIONALE ZUSAM MENARBEIT AUF DEM GEBIET DES 
PATENTWESENS (PCT) VEROFFENTLICHTE INTERNATIONALE ANMELDUNG 



(19) Weltorganisation fur geistiges Eigentum 
Internationales Biiro 

(43) Internationales Veroffentlichungsdatum 
25. April 2002 (25.04.2002) 




PCT 



II 



(10) Internationale Verdffentlichungsnummer 

WO 02/33455 A2 



(51) Internationale Patentktassifikation 7 : G02B 6/00 

(21) Internationales Aktenzeichen: PCT/DE0 1/03962 

(22) Internationales Anmeldedatum: 

17. Oktober2001 (17.10.2001) 



(25) Einreichungssprache: 

(26) VerCffentlichungssprache: 



Deutsch 
Deutsch 



(30) Angaben zur Prior itat: 

100 51 405.7 17. Oktober 2000 (17.10.2000) DE 

(71) An m elder (fur alle Bestimmungsstaaten mil A usnahme von 
US): SIEMENS AKTIENGESELLSCHAFT [DE/DE]; 
Wittelsbacherplatz 2, 80333 Miinchen (DE). 

(72) Erfinder; und 

(75) Erfinder/Anmelder (nurfur US): HIMMLER, Andreas 



[DE/DE]; Pohlweg 34, 33098 Paderbom (DE). BIER- 
HOFF, Thomas [DE/DE]; Liitersheimer Str. 22, 34471 
Volkmarsen (DE). 

(74) Gemeinsamer Vertreter: SIEMENS AKTIENGE- 
SELLSCHAFT; Postfach 22 16 34, 80506 Miinchen 
(DE). 

(81) Bestimmungsstaaten (national): JP, US. 

(84) Best im m u ngsstaaten (regional) : europaisches Patent (AT, 
BE, CH, CY, DE, DK, ES, FI, FR, GB, GR, IE, IT, LU, MC, 
NL, PT, SE, TR). 

Erklarungen gemaB Regel 4.17: 

— hinsichtlich der Berechtigung des Anmelders, ein Patent zu 
beantragen und zu erhalten (Regel 4.17 Ziffer ii) fur die fol~ 
genden Bestimmungsstaaten JP, europaisches Patent (AT, 
BE, CH, CY, DE, DK, ES, FJ, FR, GB f GR, IE, IT, LU, MC, 
NL, PT, SE, TR) 

[Fortsetzung auf der nachsten SeiteJ 



(54) Title: RAY TRACING IN HIGHLY MULTIMODAL CHANNEL WAVEGUIDES 

(54) Bezeichnung: STRAHLVERFOLGUNG IN HOCH MULTTMODALEN KANALWELLENLEITERN 



23a 34 



< 



o 




(57) Abstract: The invention relates to a method for computationally determining optical properties of a channel waveguide, accord- 
ing to which, for an incident ray entering an entrance surface, the distribution of intensity over an emergence surface is determined 
by adjoined partial rays that possibly split up. From an incident partial ray striking the outer surface, a predetermined algorithm 
determines a reflected main ray and, in so far as it is necessary according to the ensuing course of events, determines a number of 
scattered rays of a higher order, which depict partial rays and are recursively traced further. Each partial ray is either a primary ray or 
a secondary ray. When the incident partial ray is a primary ray, the reflected main ray is a primary ray once again, and the scattered 
rays are secondary rays. During the reflection of a secondary ray, only the reflected main ray is considered as the secondary ray. 
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(57) Zusammenfassung: Verfahren zur rechnerischen Bestimmung optischer Eigenschaften eines Kanalwellenleiters, bei dem fur 
einen in eine Eintrittsflache einfallenden Strahl die Intensitatsverteilung iiber eine Austrittsflache durch aneinander anschlie&ende, 
sich moglicherweise aufteilende Teilstrahlen bestimmt wird. Aus einem auf die Mantelflache fallenden Teilsu^hl wird durch ei- 
nen vorgegebenen Algorithmus ein reflektierter Hauptstrahl und, sofern nach dem folgenden Ablauf erforderlich, eine Anzahl von 
Streustrahlen hoherer Ordnung bestimmt, die ihrerseits Teilstrahlen darstellen und rekursiv weiterverfolgt werden. Jeder Teilstrahl 
ist entweder ein Primarstrahl oder ein Sekundarstrahl. Wenn der einfallende Teilstrahl ein Primarstrahl ist, ist der reflektierte Haupt- 
strahl wieder ein Primarstrahl und die Streustrahlen sind Sekundarstrahlen. Bei der Reflexion eines Sekundarstrahls wird nur der 
reflektierte Hauptstrahl als Sekundarstrahl berucksichtigt. 
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Beschreibung 

Strahlverf olgung in hoch multimodalen Kanalwellenleitem 

5 Die Erfindung betrifft die. Strahlverf olgung in hoch multimo- 
dalen Kanalwellenleitem. 

Stand rl^r TRrhn-iTc 

Zur Berechnung der Strahlausbreitung in dielektrischen Kanal- 
wellenleitem, insbesondere Lichtleitern, stehen bislang wel- 

10 lenoptische Anal ys erne thoden wie die Methode der finiten Ele- 
ment e (FEM) Oder die f Beam Propagation Method 1 (BPM) zur Ver- 
fiigung. Diese sind jedoch nur dann effizient einsetzbar, wenn 
nur eine oder wenige Moden zu berucksichtigen sind und der 
Querschnitt der Wellenleiter, bezogen auf die optische Wel- 

15 lenlange, nicht allzu grofi ist. 

Fur multimodale Stuf enindex-Wellenleiter, bei denen der Quer- 
schnitt wesentlich grofier als die Wellenlange der verwendeten 
Strahlung ist, sind hingegen Strahlverf olgungen auf der Basis 
20 geometrischer Optik effizient moglich. 

Hierbei wird (in der Simulation) ein einzelner Strahl vorge- 
gebener Richtung und Polarisation in den Wellenleiter einge- 
koppelt. Dieser tritt entweder direkt am Ende des Wellenlei- 
25 ters aus oder wird an der Wand des optischen Kanals, d.h. der 
Grenzflache des Indexsprungs, gebrochen, 

Dabei teilt sich die eintreffende Strahlung in Komponenten 
auf: einen ref lektierten Hauptstrahl, einen transmittierten 
30 Hauptstrahl, eine Anzahl von ref lektierten Streustrahlen und 
eine Anzahl von transmittierten Streustrahlen. Die transmit- 
tierten Anteile sind fur die weitere Strahlverf olgung ohne 
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Belang; lediglich ihr Energieanteil geht den ref lektierten 
Strahlen verloren. 

Eine einfache Simulation berucksichtigt lediglich den reflek- 
5 tierten Hauptstrahl, d*h. die Reflexion nullter Ordnung, und 
verfolgt dessen weitere Reflexion bis zur Austrittsf lache. 
Damit kann durch individuelle Verfolgung einer grofieren An- 
zahl von Strahlen, die der Charakteristik des Senders ent- 
sprechen, das StrahlenbUndel am Ausgang des Wellenleiters be- 
10 stimmt werden. Diese Methode ist ausreichend, wenn die Refle- 
xionen nahezu ideal sind, weil die Wandung sehr glatt ist. 

Ist jedoch die Wandung nicht glatt, so fuhrt die genannte Me- 
thode zu Ergebnissen, die entsprechenden Messungen nur 
15 schlecht entsprechen. Die Berticksichtigung der ref lektierten 
Streustrahlung jedoch fuhrt zu einem Rechenaufwand, der expo- 
nentiell mit der Anzahl der Reflexionen steigt. 

Aufgabe der Erfindung ist es daher, eine Methode anzugeben, 
20 bei der auch die Streustrahlung berucksichtigt wird, dessen 
Rechenaufwand aber allenfalls quadratisch mit der Anzahl der 
Reflexionen steigt ♦ 

Die Erfindung lost diese Aufgabe, indem sie zwischen Primar- 
25 und Sekundar strahlen unterscheidet und bei Sekundarstrahlen 
nur noch Reflexionen nullter Ordnung berucksichtigt, 

Es handelt sich um ein Verfahren zur rechnerischen Bestimmung 
optischer Eigenschaf ten eines Kanalwellenleiters, bei dem fur 

30 einen in eine Eintrittsf lache einfallenden Strahl die Inten- 
sitatsverteilung iiber eine Austrittsf lache durch aneinander 
anschliefiende, sich moglicherweise aufteilende Teilstrahlen 
bestimmt wird- Aus einem auf die Mantelflache fallenden Teil- 
strahl wird durch einen vorgegebenen Algorithmus eine reflek- 

35 tierter Hauptstrahl und, sofern nach dem folgenden Ablauf er- 
forderlich, eine Anzahl von Streustrahlen hoherer Ordnung be- 
stimmt, die ihrerseits Teilstrahlen darstellen und rekursiv 
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weiterverfolgt werden. Jeder Teilstrahl ist entweder ein Pri- 
marstrahl oder ein Sekundarstrahl. Wenn der einfallende Teil- 
strahl ein Primarstrahl ist, ist der reflektierte Hauptstrahl 
wieder ein Primarstrahl und die Streustrahlen sind Sekundar- 
5 strahlen. Bei der Reflexion eines Sekundarstrahls wird nur 
der reflektierte Hauptstrahl als Sekundarstrahl berUcksich- 
tigt. 

Kui^besidxc ei bung dar 7,<*ichnnncren 
Es zeigen 

10 Fig- 1 eine zweidimensionale Darstellung eines Modells 
der Strahlverf olgung bei einem Kanalwellenleiter 
mit rauher Oberflache, 

Fig. 2 eine zweidimensionale Darstellung, reduziert auf 
die in der Erfindung verwendeten Teilstrahlen. 



Fig. 1 zeigt eine zweidimensionale Darstellung zur Erlaute- 
rung des bisher verwendeten Verfahrens, wie es in dem Artikel 
"An Efficient Monte Carlo Based Ray Tracing Technique for the 
Characterization of Highly Multimode Dielectric Waveguides 
20 with Rough Surfaces" von Th. Bierhoff et al . , Proc. 30th Eu- 
ropean Microwave Conference, Vol.1 pp 378-381, Paris, 3rd-5th 
October 2 000, beschrieben ist. 

In Fig. 1 ist ein Langsschnitt durch einen Kanalwellenleiter 
25 20 skizziert, der einen Kern 21 und eine Umhiillung 22, auch 
als 'cladding 1 bezeichnet, hat. Der (nicht gezeigte) Quer- 
schnitt kann rund oder rechteckig sein. Andere Querschnitte 
sind gleichfalls verwendbar, aber weniger praktisch in der 
Her stel lung. Es handelt sich um e.inen Stufenindex- 
30 Wellenleiter, weil ein relativ scharfer Ubergang zwischen dem 
Kern 21 und der Umhiillung 22 vorhanden ist, so dafi eine Stufe 
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im Br echungs index entsteht, wobei bekanntermafien der Bre- 
chungsindex des Kerns 21 grofier ist als der der Umhullung 22. 
Dadurch entstehen Grenzflachen 23a und 23b, an denen ein 
Strahl im Kern 21, der auf die Grenzflache trifft, reflek- 
tiert wird. Die Grenzflachen sind bier gewellt dargestellt, 
um anzudeuten, dafi eine Oberf lachenrauhigkeit von Bedeutung 
ist; in dem behandelten Verfahren werden sie jedoch zunachst 
als Ebenen oder sonstige einfach, insbesondere analytisch zu 
beschreibende Flachen angenommen . 



10 



Der Kanalwellenleiter hat eine Eintrittsf lache 24 und eine 
Austrittsf lache 25. In die Eintrittsf lache 24 fallt ein 
(Licht-) Strahl 2 6, der an der Eintrittsf lache 24 gebrochen 
wird und als Eingangsstrahl 27 fortgesetzt wird. Ob und wie 
15 groB die Brechung ist, hangt von der Art und dem Medium der 
Ankopplung ab; fur das weitere Vorgehen ist lediglich der 
Eingangsstrahl 27 von Bedeutung. Dieser erreicht die Grenz- 
flache 23b an einem ersten Ref lexionspunkt 30 und wird dort 
gebrochen und reflektiert. 



20 



Die Berechnung eines solchen Strahlengangs mit den Mitteln 
der geometrischen Optik ist weithin bekannt. Hat der Kanal- 
wellenleiter rechteckigen Querschnitt, dann sind die Grenz- 
flachen 23a, 23b wiederum rechteckig, so dafi bei gegebener 

25 Gerade fur den Eingangsstrahl . 27 der Schnittpunkt der Geraden 
mit einer der Grenzflachen, hier der Grenzflache 23b, mit den 
bekannten Methoden der analytischen Geometrie bestimmt werden 
kann. Ist der Querschnitt rund, so ist die Grenzflache 23a, 
23b ein Zylinder, zu dem der Schnittpunkt auch noch einfach 

30 berechnet werden kann. Gegebenenf alls muB bei ungewdhnlichen 
Oder variablen Querschnitten auf iterative numerische Metho- 
den zurtickgegrif f en werden. 

Die Fortsetzung des am ersten Ref lexionspunkt 30 eintreffen- 
35 den Eingangsstrahls 30 richtet sich nach den Parametern der 
Simulation. Danach wird die Intensitat und die Richtung des 
direkt ref lektierten Strahls 31a bestimmt. Die Richtung rich- 
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tet sich nach der Tangente der idealisierten Grenzflache 23b/ 
die als rauh angesehen wird. Ferner werden diffuse Streu- 
strahlen bestimmt, die in Fig. 2 durch einen einzigen Streu- 
strahl 31b reprasentiert werden. Hierbei wird ein Raumwinkel, 
5 in dem Streustrahlung kontinuierlich vorhanden ist, zu einem 
einzigen Strahl in der Modellierung zusammengef aBt . Der 
Streustrahl 31b reprasentiert damit die Streustrahlung in ei- 
nem gegebenen Raumwinkel/ die durch eine ebene Welle mit de- 
finierter Polarisation transportiert wird; fur benachbarte 

10 Raumwinkel werden entsprechend weitere Streustrahlen be- 

stimmt, die auf gleiche Art weiter behandelt werden, aber der 
Obersichtlichkeit halber nicht dargestellt sind. Gleichfalls 
angedeutet durch Pfeile sind ein gebrochener Primarstrahl und 
ein gebrochener Streustrahl. Deren weiterer Weg ist ohne Be- 

15 deutung; lediglich ist f estzuhalten, dafi die Summe aller In- 
tensitaten gleich der des eintref f enden Eingangsstrahls 27 
sein mufi. 

Der direkt reflektierte Hauptstrahl 31a trifft an der Stelle 
20 32 wieder auf eine Grenzflache, hier die Grenzflache 23a, und 
wird dort zu einem ref lektierten Hauptstrahl 33a und mehreren 
ref lektierten Nebenstrahlen 33b, von denen wieder nur einer 
dargestellt ist. Die transmittierten Strahlen werden wieder 
nur ihrem Anteil nach beriicksichtigt . 

25 

Entsprechend trifft der reflektierte Streustrahl an der Stel- 
le 36 auf die Grenzflache 23b und ergibt einen Hauptstrahl 
und mehrere Streustrahlen. 

30 Ersichtlich ist dieses Verhalten von der Struktur durch re- 
kursive Verfahren relativ einfach programmierbar . Es ist je- 
doch auch sofort ersichtlich, daB insbesondere bei mehreren 
Streustrahlen die Anzahl der Strahlen exponentiell steigt. 
Dabei werden Strahlen, deren Intensitat unter einer vorgege- 

35 benen Schwelle liegen, naturlich nicht weiter beriicksichtigt . 
Man mufi jedoch feststellen, daB entweder die Schwelle so hoch 
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liegt, dafi die Ergebnisse nicht gut genug sind, Oder zuviel 
Rechenzeit benotigt wird. 

In dem genannten Artikel wird daher vorgeschlagen, die Kom- 
5 plexitat durch Anwendung der Monte-Carlo-Methode zu reduzie- 
ren. Diese allgemein bekannte Methode wahlt durch einen Zu- 
fallsmechanismus einen der Streustrahlen aus und ignoriert 
die anderen. Es ergab sich zwar eine Reduktion der Rechen- 
zeit, aber insgesamt noch keine einfache, schnelle und prob- 
10 lemlos benutzbare Losung. 

Durch die Erfindung kann auf die Monte-Carlo-Methode verzich- 
tet werden und ein wesentlich besseres Ergebnis in Relation 
zum Rechenaufwand erreicht werden. Dabei werden alle nach dem 

15 Modell der Reflexion entstehenden Streustrahlen weiterver- 
folgt und nicht nur einige zufallig ausgewahlte. FUr diese 
Streustrahlen jedoch wird bei nachf olgenden Reflexionen nur 
noch der Hauptstrahl berticksichtigt . Dabei wird jedesmal das 
Modell der Reflexion berechnet, urn die Intensitat des Haupt- 

20 strahls zu ermitteln; es wird also nicht unbedingt ein einfa- 
ches Spiegelmodell fur diese Nebenstrahlen eingesetzt, obwohl 
dies auch mdglich ist. 

Dies wird erreicht, indent zu jedem zu verf olgenden Strahl zu- 
25 satzlich ein Generationsindex mitgefuhrt wird, der angibt, ob 
es sich urn einen Primarstrahl oder einen Sekundarstrahl han- 
delt. Bevorzugt wird hierzu eine ganze Zahl benutzt, die Null 
fiir einen Primarstrahl und grofier Null fur einen Sekundar- 
strahl ist. 

30 

Dies ist in Fig. 2 skizziert. Die transmittierten Anteile 
wurden ohnehin nicht mehr dargestellt. Aber auch die von dem 
sekundaren Streustrahl 31b bewirkten sekundaren Streustrahlen 
werden gar nicht mehr weiterverf olgt . Hingegen wird der se- 
35 kundare Teilstrahl 33b als von einem primaren Teilstrahl 31a 
bei der Reflexion 32 ausgehend weiterhin berucksichtigt . 
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Bei der einfachen, bislang beschriebenen Ausfuhrung wird le- 
diglich ein binarer Generationsindex benutzt bzw. nur zwi- 
schen Primar- und Sekundarstrahlen unterschieden. Es ist je- 
doch auch moglich, die Simulation zu verbessern, indent der 
5 Generationsindex als ganze Zahl ausgenutzt wird. Dabei bleibt 
der Generationsindex bei dem Hauptstrahl einer Reflexion, 
d.h. der Reflexion nullter Ordnung, gleich. Bei Streustrah- 
len, d.h. Reflexionen hoherer Ordnung, wird die Ordnungsnum- 
mer zu dem Generationsindex addiert, Ob ein Strahl weiterver- 
10 folgt wird, wird dann tiber einen Schwellwert ftir den Genera- 
tionsindex entschieden. In dem oben genannten einfachen Bei- 
spiel werden uberhaupt nur Reflexionen erster Ordnung beriick- 
sichtigt, und der Schwellwert ist Eins. 

15 Durch den Gewinn an Effizienz konnte bei der Reflexion auch 
die Polarisation berticksichtigt werden, was bei der aufge- 
fuhrten Monte-Carlo-Methode ignoriert werden mufite. 



BNSDOCID: <WO 0233455A2_I_> 



WO 02/33455 

8 



PCT/DE01/03962 



Patentanspruche 

l.Verfahren zur rechnerischen Bestimmung optischer Eigen- 
schaften eines Kanalwellenleiters (20) mit einer Mantel- 
f lache, wobei fur einen in eine Eintrittsf lache (24) ein- 
5 fallenden Strahl (27) die resultierende Strahlverteilung 

uber eine Austrittsf lache (25) bestimmt wird, mit den 
Merkmalen: 

- zu einem auf die Mantelflache (23a, 23b) fallenden 
Teilstrahl (31a, 31b, 33a, 33b, 35a, 35b) wird ein re- 

10 flektierter Hauptstrahl (31a, 33a, 35a) der Ordnungs- 

zahl Null und, sofern nach dem folgenden Ablauf erf or- 
derlich, eine Anzahl von Streustrahlen (31b, 33b, 35b) 
mit einer Ordnungszahl groSer Null bestimmt, die ih- 
rerseits Teilstrahlen darstellen und rekursiv weiter- 

15 verfolgt werden, wobei der einfallende Strahl (27) der 

erste Teilstrahl ist, 
dadurch gekennzeichnet , dafe 

- jedem Teilstrahl ein Generations index zugeordnet ist, 

- der einfallende Strahl den Generations index Null hat, 
20 - bei der Reflexion der Generations index eines reflek- 

tierten Teilstrahls gleich dem um die Ordnungszahl der 
Reflexion erhohten Generations index des einfallenden 
Teilstrahls ist, 

- nur solche Teilstrahlen weiterverf olgt werden, deren 
25 Generations index eine vorgegebene Schranke nicht iiber- 

schreitet . 

2. Verfahren nach Anspruch 1, wobei die vorgegebene Schranke 
gleich Eins ist und bei der Berechnung der Reflexion al- 
lenfalls nur Streustrahlen der Ordnungszahl Eins anfal- 
30 len. 
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(57) Abstract: The invention relates to a method for computationally determining optical properties of a channel waveguide, accord- 
ing to which, for an incident ray entering an entrance surface, the distribution of intensity over an emergence surface is determined 
by adjoined partial rays that possibly split up. From an incident partial ray striking the outer surface, a predetermined algorithm 
determines a reflected main ray and, in so far as it is necessary according to the ensuing course of events, determines a number of 
scattered rays of a higher order, which depict partial rays and are recursively traced further. Each partial ray is either a primary ray or 
a secondary ray. When the incident partial ray is a primary ray, the reflected main ray is a primary ray once again, and the scattered 
rays are secondary rays. During the reflection of a secondary ray, only the reflected main ray is considered as the secondary ray. 
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(57) Zusammenfassung: Verfahren zur rechnerischen Bestimmung optischer Eigenschaften eines Kanalwellenleiters, be, dem fur 
einen in eine Eintrittsflache einfallenden Stnml die Intensitatsverteilung iiber eine Austritteflache durch anemander "^*"^> 
sich moglicherweise aufteilende Teilstrahlen bestimmt wird. Aus einem auf die Mantelflache fallenden Te, ^^J^* 
nen vorgegebenen Algorithmus ein reflektierter Hauptstrahl und, sofem nach dem folgenden Ablauf erf orderbch, eine Anzah von 
Streustrahlen hoherer Ordnung bestimmt, die ihrerseits Teilstrahlen darstellen und rekursiv weiterverfolgt werden Jeder Tertstrahl 
fstenwedereinPrima^ 

strahl wieder ein Primarstrahl und die Streustrahlen sind Sekundarstrahlen. Bei der Reflexion eines Sekundarstrahls w,rd nur der 
reflektierte Hauptstrahl als Sekundarstrahl beriicksichtigt. 
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